Homologous Expansion and Contraction 


Many problems in stellar structure can be addressed by approximat- 
ing the expansion (or contraction) of the star as homologous. This is 
equivalent to saying that the fractional rate of change in the star is a 
constant, t.e., 
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(An example of such a system would be an expanding polytrope, which 
keeps the same value of n. Since M œx € x r/rn, the mass elements 
of the star will remain at homologous points; only the scale length rn 
changes.) 


Although homologous expansion (or contraction) does not often occur 
(after all, it does require that the whole star act in unison), it is a good 
first-order approximation for many problems, and is useful for quick 
analyses. For instance, we can estimate the energy released (or ab- 
sorbed) when a star homologously contracts (or expands). According 


to (20.1.1) 
d dlnr 


If we exchange the two derivatives, then 
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Setting this to zero then implies that for homologous movement 
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Similarly, we can easily compute how the hydrostatic pressure changes 
for an homologous mass point under expansion. From (2.2.4) 
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which for constant t/r implies 
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Combining (20.1.3) and (20.1.4) with the equation of state for a chem- 
ically homogeneous gas 
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gives an expression for temperature 
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Recall now that the energy associated with the gravitational expan- 
sion/contraction of a star is 
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which, through (1.19) is 
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Substituting using the equation of homology then gives an estimate for 
the total amount of energy associated with an expansion or contraction. 
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For an ideal gas with a = 6 = 1 and Vaa = 2/5, 
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Note that the value is negative; contraction by the star releases energy, 
while expansion absorbs energy. 


The Gravothermal Specific Heat 


Another use of dynamic homology involves the calculation of 
Gravothermal Specific Heat. Two familiar thermodynamic con- 
cepts are the specific heat for a system kept at constant pressure 
(cp) and the specific heat for a system in a constant volume 
(cy). However, for some stellar applications, a more interest- 
ing quantity is the specific heat of a system kept in hydrostatic 
equilibrium. This is called the Gravothermal Specific Heat. 


To compute the gravothermal specific heat, c*, let us pick a re- 
gion of the star with radius, r, and assume that, at least locally, 
any expansion or contraction near that region is homologous. 
(That is, dr = r dz, and after expansion, the nearby shells have 
radii r + dr = r(1 + dz), where dz is the same for all shells.) 
From (20.1.3) and (20.1.4) 
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Hence, if the expansion is homologous 
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Now let’s substitute this into the equation of state 
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If the chemical composition is constant, then 
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which implies 
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Now let’s use this result in the first law of thermodynamics. 
Recall that when written in terms of cp and Vaq, the first law 


was 
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(1.21) 
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For homologous expansion (or contraction), we then have 
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Thus, the gravothermal specific heat acts like any other specific 
heat: if an amount of heat, dq is added to the system, then the 
temperature of the system will rise by dT = dq/ œ. 


Equation (20.2.6) has interesting implications for normal stars. 
For an ideal gas equation of state (a = 1,6 = 1, Vaq = 2/5), c* 
is negative. Thus, if the heat content of a core is perturbed by 
an amount dq, the temperature of the core immediately drops, 
thereby reducing the production of energy. (If this did not 
happen, then we would have a thermal runaway!) 


The reason stars are stable is that any extra energy that is 
produced goes into performing PdV work on the surroundings. 
If heat is added to the center of an (ideal gas) star, dT < 0 
(by 20.2.5), which implies dP < 0 (from 20.2.4), which implies 
dx > 0 (from 20.2.2). Thus the star expands! The extra heat 
is used to expand the star, and more energy is used in the 


expansion than was originally added. The star must extract 
the requisite additional energy from its internal energy. 


Note also that under homologous expansion or contraction, 
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Therefore, since dx is a constant, all these ratios are also con- 
stants, with 
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Stability in Core Burning 


Now let’s consider how a star will react to changes in its rate of 
nuclear reactions. In equilibrium, the energy generated inside a 
sphere with mass, M, is balanced by the luminosity transported 
outward, t.e., 
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Now, let’s perturb the energy generation by an amount de. In 
general, the dynamical timescales of stars are short, and the ther- 
mal timescales relatively long, so that during the perturbation, the 
star will remain in hydrostatic equilibrium, but be out of thermal 
balance. Thus 
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or, if we divide through by £ = eM, 
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Now let us derive an expression for d£/L in the case of radiative 
energy transfer. From (3.1.5) 
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or 
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This can be simplified. First, recall that from the definition of 
homologous movement and from (20.2.7) 
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Next, consider that the opacity, x, is itself a function of density 
and temperature, with k x pPT4. Thus 
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where we have used (20.2.8) in substituting for dp/p. Moreover, 
under homology dT’/T is a constant, so 
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Putting this all together, we get 
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or, still more simply, 
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Similarly, we can substitute for de/e using the energy generation 
coefficients. If we let € x pòT” and again use (20.2.8) then 
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If we now substitute this and (20.3.4) into (20.3.2), we get 


ô dT M „dT 
i dT K dT 
= 20.3.7 
dt eT ( ) 
where 
EA E dh (3A +3p +4) (20.3.8) 
-M i da —3 b a 


Note the implications of this equation. In the case of an ideal 
gas (œa = 6 = 1), proton-proton burning (v ~ 5, A = 1), and 
Kramers opacity (q = —7/2, p = 1), K = 7.5£/M. Since c* 
is less than zero, a positive temperature perturbation dT/T will 
result in a negative dT/dt, i.e., a net cooling. As a result, the 
temperature perturbation will be damped away. The coefficients 


for CNO burning, electron scattering, and H~ opacity make K 
even more positive, increasing the stability. 


Now consider a non-relativistic degenerate equation of state with 
helium fusion. From (7.3.7) P œ př/3, which implies 6 = 0, and 
a = 3/5. Thus, we have a positive gravothermal specific heat 
(c* = cp), and a positive value of K (v +q — 4 ~ 32). Under these 
circumstances, a positive heat fluctuation leads to an increase in 
temperature, and a further increase in the nuclear reaction rate. 
We have a thermonuclear runaway! Moreover, from (20.2.8), 
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Consequently, a stellar core can experience a thermonuclear run- 
away, while keeping the density (and, by implication, the pressure) 
constant. 


Stability in Shell Burning 


The analysis of the thermal stability of shell burning stars follows 
that for core-burning stars, with one exception. With core burning, 


the mass within a burning region is M = srr p, and so under 
homology 
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However, consider the case of a shell of thickness D at a radius ro 
which is defined by the inert core below. In this case 
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Since ro does not change, and the mass in the shell does not change, 
dD = dr, and 


M 1 — MadaD_ dD 
Agr? D? PD 


or 


E Le me N (20.3.9) 


Since the equations of state and of hydrostatic equilibrium do not 
change, we can now substitute this new expression for dp and re- 
calculate c*. The result is 
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(Obviously, this is identical to the original expression, except that 
the constant “3” has been replaced by “r/D” in the denomintor.) 


Now let’s again consider the equation of thermal stability (20.3.7). 
Since K is positive, a negative value of c* implies thermal stability, 


while a positive value denotes a thermal runaway. As demonstrated 
above, if the gas conditions are degenerate, ô = 0, and again we 
have a runaway. However, suppose we are dealing with an ideal gas 
(which is the relevant condition, since the gas in the shell is outside 
the core). For core burning stars, an ideal gas equation of state 
always produces a negative value of c*, so the star is stable. But 
according to (20.3.11), if D is small, r/D will be large enough to 
make c* positive. ‘Thus, nuclear burning in thin shells is susceptible 
to thermonuclear runaways. 


In essence what is happening is that if D/r is small, a small change 
in D causes a large change in p, but, in absolute terms, a very minor 
change in dr/r. Thus the layers above, which are in hydrostatic 
equilibrium, are barely affected, and dP/P = 0. The effective 
equation of state is therefore p ~ 1/T, for which a = 0 and 6 = 1. 
Thus, fusion under these conditions can be unstable and runaways, 
or thermal pulses, can occur. 


